Distinguishing d-wave from highly anisotropic s-wave superconductors 
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Systematic impurity doping in the Cu-0 plane of the hole-doped cuprate superconductors may 
allow one to decide between unconventional ("d-wave") and anisotropic conventional ("s-wave") 
states as possible candidates for the order parameter in these materials. We show that potential 
scattering of any strength always increases the gap minima of such s-wave states, leading to acti- 
vated behavior in temperature with characteristic impurity concentration dependence in observable 
quantities such as the penetration depth. A magnetic component to the scattering may destroy the 
energy gap and give rise to conventional gapless behavior, or lead to a nonmonotonic dependence 
of the gap on impurity concentration. We discuss how experiments constrain this analysis. 
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Introduction. A number of recent experiments on hole- 
doped cuprate superconductors have provided evidence 
for a superconducting state with very large anisotropy, 
consistent with actual gap nodes on the Fermi surface, jll- 
0] The set of experimental results indicating the existence 
of low-energy quasiparticle excitations have been inter- 
preted in terms of an unconventional, "d- wave" pairing 
state, where we use the term unconventional to mean that 
the superconducting order parameter breaks additional 
symmetries of the normal state beyond the usual gauge 
symmetry. Such an order parameter iS.k^ has a non- 
trivial phase variation over the Fermi surface and changes 
sign at the node. Since the quantities measured in these 
experiments usually depend on the order parameter only 
through the quasiparticle energy Ek = (Cfe^ + 
where is the single-particle energy measured relative 
to the Fermi level, it is easy to see that an identical result 
would be obtained by a measurement on a hypothetical 
state with order parameter A^" — \Ak'^\, which would 
vanish at the same nodal points but never change sign. 
Since the nodal points in this case are accidental rather 
than being enforced by symmetry, it is more realistic to 
consider a highly anisotropic s-wave state with very deep 
gap minima but no nodes. Such an order parameter has 
in fact been proposed by Chakravarty et al.,[l| and is 
quite difficult to distinguish from a similar d-wave state 
if the experiment does not measure temperatures sub- 
stantially below the gap minima and is not sensitive to 
the gap phase variation. 

As a consequence of these ambiguities, methods 
of distinguishing between unconventional states and 
anisotropic conventional states are of great importance. 
Josephson tunneling experiments are sensitive to the or- 
der parameter phase and therefore in principle capable 
of deciding this question. 0, [1] At present, however, dif- 
ferent Josephson experiments of slightly different con- 
cept and design have reached differing conclusions re- 
garding the order parameter symmetry[3,H, We have 
therefore reexamined the well-studied problem of dirty 
superconductors with an eye towards designing further 
tests which may be capable of distinguishing d-wave and 
highly anisotropic s-wave states. We find that systematic 
impurity doping experiments are indeed sensitive to the 
order parameter phase, albeit indirectly, and can provide 
important evidence towards the resolution of this ques- 



tion. 

Model — Potential Scattering. For illustration's sake 
we consider a dr^2_y2 state over a cylindrical Fermi sur- 
face, Afe = Ao cos 2(j} and a hypothetical s-wave state 
|AoCOs20|. Norman |14| has shown that weak poten- 
tial scatterers eliminate the nodes in the s-wave case 
at the points (j) = 7r/4, 37r/4, increasing the gap in 
these directions monotonicall y w ith impurity concentra- 
tion. It was in fact argued in [14'| that the dependence of 
angle-resolved photoemission (ARPES) data on sample 
aging [H may be construed as evidence for s-wave super- 
conductivity, but there are alternative explanations pe- 
culiar to the ARPES configuration. Q Here we consider 
these two states further, and investigate the effects of 
strong scattering as well as spin scattering, and try to 
make predictions for bulk thermodynamic and transport 
experiments. 

The properties of a dx2_y2 state in the presence of 
elastic impurity scattering have been extensively inves- 
tigated in recent months. |l2l. [l5| but are in fact generic 
to states with lines of nodes on the Fermi surface in 3D 
(point nodes in 2D) investigated in the context of heavy 
fermion superconductivity. [igI [17| An infinitesimal num- 
ber of impurities suffice to make the density of states at 
the Fermi level nonzero, [3| giving rise at low tempera- 
tures T <C Tc to contributions which vary with tempera- 
ture as their normal state analogues, but with a smaller 
prefactor which scales with impurity concentration. The 
penetration depth, which does not have a normal state 
analogue but varies as ^ T/Aq in the pure d-wave state, 
is known to cross over to a behavior in this so- called 
"gapless" regime. [ni[il At the same time, the actual en- 
ergy gap in the angle-resolved density of states remains 
zero along the nodal directions. All these characteristics 
may be understood as consequences of the exact vanish- 
ing of the anomalous impurity self-energy which occurs 
in most — but not all — unconventional states. 

The essential differences between s- and d-wave states 
may be understood by examining the single particle ma- 
trix propagator g averaged over impurity positions, given 
by:18| 



|A.. n (1) 



where the are the Pauli matrices and At. is assumed 
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to be a unitary order parameter of s- or d-wave type in 
particle-hole and spin space. The propagator (1) has the 
form of the propagator for the pure system with renor- 
malized frequency w = cj — I]o(w), single-particle energy 
£,k = £,k + ^3(0;), and order parameter = Afc -I- Ei(w), 
where the self-energy due to s-wave impurity scattering 
has been written E = Sir'. For the particle-hole sym- 
metric systems we consider here, renormalization of the 
single-particle energies can be important for arbitrary 
scattering strengths, but are small in either the weak 
or strong scattering limit. We will neglect them in what 
follows. 

As alluded to above, in odd-parity states and states 
with certain reflection symmetries like dx2-y2, the off- 
diagonal self-energy Ei vanishes identically and the gap 
is unrenormalized, (A^ = Aj.). Potential scatter- 
ers are then pairbreaking, in "violation" of Anderson's 
theorem, [li'l but the angular (e.g., nodal) structure of 
the gap is not changed. By contrast, in the anisotropic 
s-wave case the order parameter A^ is always renormal- 
ized by a positive shift which is independent of k in the 
s-wave scattering approximation. This leads to a smear- 
ing of the energy gap anisotropy leading eventually to 
an asymptotically isotropic gap in the dirty limit, as im- 
plied by Anderson (Toj and calculated explicitly by various 
authors. |20| 

In the absence of renormalizations, the self- 
energies are given in a t-matrix approximation by Eq — 
TGo/D; El = TGi/D, where F = nJ{nNo) is a scat- 
tering rate parameter depending only on the concentra- 
tion of defects Ui and the density of states at the Fermi 
energy, Nq, while the strength of a single scattering is 
characterized by the cotangent of the scattering phase 
shift, c. Here D = + Gi^ — Go^ is the denomi- 
nator determining the bound state spectrum, and the 
Ga = (l/27rA'o)EfcTr[r"(7] are components of the inte- 
grated, disorder-averaged propagator. The Born limit 
corresponds to c 1, so that F/c^ ~ Tfq = F/(l -|- c^), 
where Fjv is the scattering rate in the normal state due 
to impurities. The unitarity or strong scattering limit 
corresponds to c = 0. 

Order parameter, critical temperature, and energy 
gap. We first solve the Dyson equation for the renor- 
malized propagator (1) together with the gap equa- 
tion, A(fc) - rj:„^^, T4fc'Tr(Ti/2)g(fc',tj„), where 
Vkk' = Vd.s^d.s{k)^d,s{k') is the phenomenological pair 
interaction assumed. The order parameter is Aj. = 
Ao'^'*<&d,s(^), with ^d.s — cos 2(/), I cos 2(/)| for d and s 
wave, respectively. The initial slope of Tc suppression, 
dTjdVN = where x = - {<^s?]/{^sf is 

1 — 8 /tt^ for the s-wave and 1 for the d-wave state consid- 
ered. In the d-wave case the critical temperature contin- 
ues to drop rapidly to zero at a critical concentration of 
n'l = '7T'^NoTco/2e'' , whereas the decrease becomes more 
gradual as the gap is smeared out in the s-wave case, 
finally varying^ as T, - r,o[l - xln(1.154Fw/^Teo)]. 

It is important to recognize that the renormalized or- 
der parameter A^ in the s-wave case is only indirectly re- 
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FIG. 1: Normalized density of states N{u})/No for s- and d- 
wave order parameters vs. reduced frequency oi/Aq, shown for 
various potential scattering rates Pjv / Ao in Born approxima- 
tion. 
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FIG. 2: Normalized density of states N{iu)/No for and d- 
wave order parameters vs. reduced frequency lj/Aq, shown 
for various potential scattering rates F/Aq in unitarity limit, 
c = 0. 



lated to the actual energy gap ilc in the system, given by 
the maximum frequency uj such that the angle-resolved 
density of states N{ls.,uj) = Im Tr [g(k, aj)]/7r = for 
all k. A simple estimate shows that for small scattering 
rates, flc ~ r(r^ in Born limit). In the dirty limit F — 
00, the s-wave superconductor becomes isotropic with a 



BCS density of states N{uj) = Re [uj/{uj^-Aa 
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shown in Fig. 1. In contrast to a d-wave superconductor, 
the self-energies obtained in the Born approximation and 
in the resonant scattering limit are almost equivalent in 
the highly anisotropic s-wave system. This insensitivity 
to larger phase shifts arises because of off-diagonal self- 
energy corrections which prevent the occurence of poles 

in the t-matrix, - Gq^ + Gi^ ~ 0(1) for ah c<. Den- 
sities of states for both types of states in the limit of 
resonant scattering are shown in Fig. 2. 

London penetration depth. The opening of the energy 
gap with increasing impurity concentration is an indeli- 
ble signature of s-wave superconductivity. It will obvi- 
ously give rise to activated behavior for T <IC fig in a 
wide range of thermodynamic properties, of which we 
have chosen to discuss only one for purposes of illustra- 
tion, the temperature-dependent magnetic penetration 
depth. For the model states and Fermi surface under 
consideration, this may be expressed as [Ao/A(T)]^ = 
/ duj tanh/3w/2 / #/27rRe Al/{u^ - Klf'^ where Aq 
is the pure London result at T = 0. The penetration 
depth in a d-wave superconductor (Fig. 3, bottom half) 
is known to vary as A(T) ~ Aq -I- C2T^ at the lowest 
temperatures. [l ll [T2j over a temperature range which 
widens with increasing impurity concentration. The co- 
efficient C2 decreases, as F"^ in the Born limit and F~^/^ 
in the resonant scattering case. The corresponding acti- 
vated behavior in the anisotropic s-wave case is easy to 
distinguish from the d-wave case when plotted against 
(T/Tc)^ as also shown in Fig. 3. The important exper- 
imentally relevant signature is of course not simply the 
exponential behavior, but the increase in the activation 
gap with impurity concentration. 

Spin scattering. A simple defect like a vacancy or Zn 
ion in the Cu02 plane may not behave simply as a poten- 
tial scatterer, as assumed above. In the presence of large 
local Coulomb interactions, a magnetic moment may 
form around the defect site, giving rise to spin-flip scat- 
tering of conduction electrons. |13| This poses the most 
serious obstacle for the direct application of the principle 
distinguishing d-wave from anisotropic s-wave systems 
outlined above, since magnetic scattering will lead to gap- 
less superconductivity as in the usual Abrikosov-Gor'kov 
theory. Furthermore, even if a gap remains, strong spin- 
flip scattering may lead to bound states within it[2ll. [22j 
which may give rise under the proper circumstances to 
a residual density of states N{uj 0) as in the d-wave 
case. Here we investigate the competition between the 
opening of the energy gap in the s-wave state due to 
potential scattering and gapless behavior due to mag- 
netic scattering. To this end we add a term JS ■ a to the 
Hamiltonian, where S is a classical spin representing the 
impurity and a is the conduction electron spin density, 
and study the system in an average t-matrix approxi- 
mation analogous to the one applied to the pure poten- 
tial scattering case. The self-energies found in the pres- 
ence of both types of scattering reduce in the isotropic 
s-wave case to those given by Shiba,[22| but are com- 



\ 
O 
r< 



O 



0.8 



0.6 



0.4 



0.8 



0.6 



0.4 



\r/Tco=D 


S 


\ 


c=0 


1 




2 




; ' 

\ 


d 




c=0 



















0.1 0.2 

FIG. 3: Temperature dependence of normalized magnetic 
penetration depth (Ao/A(r))^ for s- and d- wave order pa- 
rameters vs. reduced temperature (T/Tc)^, shown for various 
potential scattering rates V /Tco in unitarity limit, c = Q. 
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FIG. 4; Induced energy gap normalized to clean gap maxi- 
mum, fic/Ao, vs. potential scattering rate Tm/Tco for differ- 
ent ratios, V%/Vn of magnetic to potential scattering rates. 



plicated and will be given elsewhere. We find that un- 
til the dimensionless exchange JNq becomes of 
the results for the s-wave system are very similar to 
those obtained in the simpler Born approximation, as 
discussed above. In this case, Eq = (Fat + F|^)Go and 
Si = (_rjv-HF^)Gi,[ll where F^ = n,J^ S{S + 1)t:No. 
The induced gap, VLq, in the s-wave system may then be 
shown to vary as VLq ~ Fat — T^^ > 0, but the effects of 
self-consistency rapidly become important as the concen- 
tration is increased. In Fig. 4, we plot Q.q as a function 
of the impurity concentration through the parameter F^r 
for various assumptions about the scattering character 
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of the impurity ion, where the quantity /T^ specifies 
the relative amount of magnetic scattering. The destruc- 
tion of the induced gap takes place because the system 
becomes insensitive to large amounts of potential scatter- 
ing, but magnetic impurities continue to break pairs even 
at large concentrations. The gap is nevertheless found to 
persist into the very dirty limit even for systems where 
the magnetic scattering is nearly as strong as the poten- 
tial scattering. 

For weak spin scattering, the bound state in the t- 
matrix approximation is found to lie at w 3> just 
below the average order parameter Aavg deep in the con- 
tinuum, and thus plays no role. Stronger spin scattering 
does not change this qualitative behavior at low concen- 
trations until JNq ~ 1 when the bound state lies at the 
Fermi level in the classical spin approximation. |22j In this 
case the Kondo effect, neglected here, also becomes im- 
portant. It is known from other analyses ^235 that the 
bound state lies near the Fermi level, and will therefore 
give rise to a residual density of states N{uj — > 0), only 
when Tk — Tc- For any other ratio of Tk/Tc, the bound 
state will lie at an energy corresponding to an apprecia- 
ble fraction of the average gap in the system, and hence 
be irrelevant for our purposes. 

Clearly a quantitative estimate of the relative size of 
r and F^ is required to decide whether spin scattering 
plays a role in real high-Tc materials with simple defects. 
Walstedt and co-workers estimated JNq ~ 0.015 for a Zn 
ion in YBCO, implying that Zn is a nearly pure potential 
scatterer in this system. [23| On the other hand, Maliajan 
et al.jlll estimate JNq ~ 0.45. For a 1% Zn concentra- 
tion, a magnetic moment of 0.36 /is for Zn in fully oxy- 
genated YBCO[25] and a density of states of 1.5/eyj25|. 
we find F^ ~ 1 x 10^"* eV. From the residual resistivities 
of Zn-doped YBCO crystals, [1^ we estimate that a 1% 
Zn sample corresponds to a total impurity scattering rate 
of F™^' ~ 1 X 2 X 10 '^eV. assuming that the inelastic 
and elastic contributions to the scattering rate add in- 
coherently. This suggests that potential scattering must 
dominate the total elastic rate, F|^ <C F. On the other 
hand, the large value of JNq ~ 0.45 deduced for a Zn 
ion[25| means that the Kondo effect may be important, 
and that we cannot completely rule out the possibility 
that a bound state sits very close to the Fermi level. 



Conclusions. There is by now a considerable body of 
experimental data supporting the picture of gapless su- 
perconductivity in the cuprate high-Tc materials, with a 
residual density of states and low-temperature behavior 
varying qualitatively according to the d-wave plus res- 
onant scattering model, [ll, This data stands in ap- 
parent contradiction to the well-known effect of small 
amounts of potential scatterers on anisotropic s-wave 
superconductors, namely the smearing of energy gap 
anisotropy. This continues to hold even for extremely 
anisotropic systems with nodes, as illustrated by the sim- 
ple theory presented here for a representative order pa- 
rameter. We believe that this data strongly suggests that 
the pairing is unconventional in these materials, but the 
above analysis does not as it stands allow one to dis- 
tinguish among possible candidate unconventional states 
(e.g,, dx2-y2 and dxz) without further quantitative com- 
parison. It should be noted that time-reversal breaking 
unconventional states with a gap will become gapless in 
the presence of pure potential scattering. 

As we have briefly discussed, the major difficulty in- 
herent in such an analysis is the possibility that even 
an apparently "inert" impurity such as Zn or a vacancy 
in the Cu-0 planes may induce local spin correlations 
in the strongly interacting electron system, leading to 
spin-flip scattering. Ruling out gapless superconductiv- 
ity induced by magnetic scattering then becomes a quan- 
titative problem. Gapless behavior in films suggests that 
a resonant scattering mechanism of some type must be 
present in order to induce a significant residual density 
of states with comparatively little Tc suppression. We 
have shown, however, that resonant potential scattering 
does not take place in s-wave systems, and argued that 
low-energy resonant spin scattering is much less likely 
than in the isotropic case. We have furthermore made a 
crude estimate of the importance of spin-flip scattering 
in Zn-doped YBCO crystals which indicates these ma- 
terials are dominated by potential scattering and should 
therefore exhibit an induced gap if the superconducting 
state is s-wave. 

We are grateful to N. Goldenfeld for useful discussions 
and encouragement. Some numerical computations were 
performed on the Cray-YMP at the Florida State Uni- 
versity Computing Center. 
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